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When the background density in a bounded plasma is modulated in time, discrete modes be¬ 
come coupled. Interestingly, for appropriately chosen modulations, the average plasmon energy 
might be made to grow in a ladder-like manner, achieving up-conversion or down-conversion of the 
plasmon energy. This reversible process is identified as a classical analog of the effect known as 
quantum ladder climbing, so that the efficiency and the rate of this process can be written imme¬ 
diately by analogy to a quantum particle in a box. In the limit of densely spaced spectrum, ladder 
climbing transforms into continuous autoresonance; plasmons may then be manipulated by chirped 
background modulations much like electrons are autoresonantly manipulated by chirped fields. By 
formulating the wave dynamics within a universal Lagrangian framework, similar ladder climbing 
and autoresonance effects are predicted to be achievable with general linear waves in both plasma 
and other media. 

PACS numbers: 52.35.-g, 52.35.Mw, 42.65.-k, 47.10.Df 


Introduction. — Quantum mechanics is well known 
to be closely related to the mechanics of classical waves 
0-[3[. This permits applying common techniques for 
manipulating quantum and classical systems and helps 
bridging seemingly different areas of physics. One im¬ 
portant technique to study in this context is ladder 
climbing (LC), which is the successive transfer of quanta 
through nonequally-spaced energy levels due to an oscil¬ 
lating driving force with chirped frequency 00- The 
system energy changes with time in a ladder-like manner 
during LC, with each transition described by the famous 
Landau-Zener (LZ) theory 0. In the limit of continuous 
spectra, the effect has been widely known as classical 
autoresonance (AR), enjoying numerous applications in 
phy sics of plasmas 00, fluids 0, Josephson junctions 
TlOj . optics 0, and even planetary dynamics [3]. In 
contrast, the discrete nature of LC is visible only in sys¬ 
tems with sufficiently discrete spectra and, so far, has 
been studied exclusively in quantum contexts 0 i, Il3l — 
[l5| . Whether classical systems can exhibit LC has re¬ 
mained an open question. 

Here we report the first theoretical prediction of LC 
in a classical system, namely, in an ensemble of plasma 
waves. For simplicity, we consider one-dimensional col¬ 
lisionless plasma, with nondissipative Langmuir waves, 
whose spectrum is quantized due to the boundary con¬ 
ditions. We derive a Schrodinger-type equation for the 
“plasmon wave function”, which is a classical measure 
of the electric field whose norm (the total wave action) 
is manifestly conserved under mode coupling. This sys¬ 
tem is mathematically equivalent to a quantum parti¬ 
cle in electrostatic potential. Hence, plasmons can be 
manipulated by resonant modulation of the underlying 
medium, much like electrons and molecules are manipu¬ 
lated by resonant external fields [30. In particular, 
we show that plasmons can exhibit both LC and AR and 
can be controllably transported up and down in momen¬ 


tum space. Finally, we report a unifying Lagrangian for¬ 
mulation of the problem that paves the way for applying 
these techniques to general classical waves. 

Basic equations. — For simplicity, we consider an elec¬ 
tron plasma described by a hydrodynamic model. The 
equations for the electron density n e , electron flow veloc¬ 
ity u e , and the electric field E are then as follows: 
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Here, — e and m e are the electron charge and mass, P is 
the electron pressure, Ze is the ion charge, and nt is the 
ion density. We neglect high-frequency oscillations of m 
and consider m to be a slow function, Zn* = no+nd(t, x); 
here no is the unperturbed electron density, and nd is a 
prescribed driving modulation. Such a modulation can 
be created by external fields, e.g., by means of pondero- 
motive forces. Then, n e = no + rid(t,x) + n(t,x ), where 
n <C no determines a small uncompensated charge den¬ 
sity due to electron inertia. For simplicity, we adopt 
an isentropic model, P = P(n e ) ~ Po + 3 m e v th (n + 
rid) + P(n + nd) 2 / 2; here v t h is the electron thermal 
speed, and R = d 2 P/dn 2 \ no , which are considered con¬ 
stant. We assume lrard-wall boundary conditions, so 
u e |x=o = u e \x=t = 0, where I is the plasma length. We 
also assume £ , | a; =o = E\ x —i = 0, so any field is rep¬ 
resentable as a series of sin (k m x), where k m = rtm/l. 
Then, according to Eqs. O-©, the boundary conditions 
for the density must be d x n\ x —o = d x n\ x —e = 0, so n is 
a series of cos(k m x). We consider the external driving 
modulation to be the Nth standing-wave mode, 


rid(x, t) = no A cos ( kNX) cos fd{t). 


( 4 ) 
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We assume that A <Cl and u>d = <p <C u> p , where 
ojp = 47re 2 no/m e . Hence, Eqs. ©-© can be combined 
in a single equation. Specifically, let us subtract the spa¬ 
tial derivative of Eq. © from the temporal derivative 
of Eq. (JT]) ; then substitute d x E = — 47ren, which flows 
from Eq. m, and integrate over x. Assuming that n is 
small enough, we neglect terms nonlinear in E = 0(n), 
and we also neglect the slow driving force, O(E 0 A 1 ), as 
nonresonant to the rapid plasma oscillations. As we es¬ 
timate below, resonant terms of higher orders in A can 
be dropped too. To the lowest order in u>d [3; this gives 
the following linear dimensionless equation for E, 
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Here and further we measure time in units ui p 1 and 

length in units v t h/u v \ also, rid = ^d/n o> and R = 
Rno/(mv 2 h ). Next, we decompose the field into unper¬ 
turbed eigenmodes, 


P = Re^P m e lUmt y/2ji sin(k m x), (6) 

m= 1 


where E m (t ) are complex coefficients. To zeroth order in 
A, Eq. © yields the dimensionless dispersion relations 
= 1 + /3m 2 , where /3 = 37 t 2 /1 2 is analogous to the an- 
harmonicity parameter in a quantum oscillator Q. We 
also assume /3m 2 <C 1, as usual. It is convenient to intro¬ 
duce new variables ijj m via Pmipm(t) = E m e~ lulmt . Here 
p m are constants such that in the unperturbed system, 
|'*/’m | 2 are the actions of individual modes; specifically 
one finds p m = ( 87 r/u m )^ 2 [ 3 - The equations for ip m , 
obtained via Fourier-transforming Eq. m , are: 

E hm,m ’1pm '1 (7) 

m‘ 

sm(k m x)J- sm(k m >x) dx, (8) 

where T = rid — d x Rrid d x is a differential operator. 
Since h m ,m' is Hermitian, the evolution of ip m is man¬ 
ifestly unitary; i.e., the wave total action \'4 > m\ 2 is 
conserved. The vector %p = (ipi, ip 2 , ■ ■ ■) can then be 
understood as the plasmon wave function in the en¬ 
ergy representation. Next, we will assume the reso¬ 
nance condition tod ~ u} mt m+N, where the level spac¬ 
ing is LUm,m+N = u m +N - ~ /3N(m + N/2). Then, 

/3 ~ Wd 1, and therefore, in the already small coupling 
term, we must adopt ~ 1 , i.e., p m ~ \/ 87 r and neglect 
the higher order terms in pressure. This leads to 

A 

ilp m = U} m 4> m + — plpm-N + V’m+Tv) COS p d - (9) 
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FIG. 1: (color online) Langmuir wave evolution in the LC 
regime (P 2 = 10). Left - spectral representation; snapshots of 
the occupation numbers, \ipm\ 2 , vs mode numbers, m at times 
r = [—5,15, 35, 55]. Right - spatial representation; snapshots 
of the wave field, E, vs x at the same times. The inset in 
subplot (a) shows the energy growing with time in a ladder¬ 
like manner, where the transitions occur at the theoretically 
predicted times r m = mP 2 - 


quantum particle in a square potential well 01. This 
can be understood as follows: at weak spatial dispersion, 
all oj m are close to w p , so Eq. m permits a quasiopti- 
cal approximation, turning into the standard quantum 
Schrodinger equation with serving as an effective po- 
tential. (This analogy has been noted, e.g., in Ref. (21].) 

LC regime. — Classical AR was previously studied 
in the infinite square potential well as a limiting case 
(j —> oo) of the potential V 0 = x 2 i 0; but quantum LC 
was not 0 ■ F° r studying LC in our system, we consider 
the modulation © with iV = 1 and a monotonically in¬ 
creasing driving frequency, Ud = £*> 1,2 + cit, where a > 0. 
At t = 0, this will initiate resonant transitions between 
levels m = 1 and m = 2, which we denote as 1 —> 2. 
Later, transitions between higher levels, m —> m + 1, oc¬ 
cur when the resonance condition, uid = Wm,?n+i, is sa E 
isfied. Following the quantum LC theory 0, Q, we define 
slow time r = y/at, driving parameter Pi = A/(4y/a), 
and anharmonicity parameter, P 2 = fi/y/a, which plays 
the role of an effective Planck constant in the classical 
system. If P 2 1 + P\ , the system remains in the quan¬ 
tum LC regime, when only two levels are resonantly cou¬ 
pled at any given time. Then, transitions to —> m + 1 
occur at times r m = mp 2 , where each such transition can 
be described by the commonly known LZ theory [§, jij 
and thus has a probability 

Vm^m +1 = 1 - exp (-7 tP 2 /2) . (10) 


Note also that, for /3m 2 1, Eq. ([9]) is equivalent to the At large enough Pi, « 1, i.e., all quanta are 

energy representation of the Schrodinger equation for a being transferred, and the dynamics in this limit is char- 
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FIG. 2: (color online) Langmuir wave evolution in the AR 
regime (P 2 = 0.1). Left - spectral representation; snapshots 
of the occupation numbers, |'0™.| 2 , are shown vs mode num¬ 
bers, m, at times r = [0,3,6, 9], illustrating an autoresonant 
plasmon trapped in a density modulation effective potential. 
Right - spatial representation; snapshots of the wave field, 
E (solid blue lines), and the driving density modulations, rid 
(dashed red lines), are shown vs x at the same times. The 
inset in subplot (a) shows the energy growing with time. 


acterized by successive two-level LZ transitions. 

In Fig. m we illustrate this LC dynamics of Langmuir 
modes by numerically simulating Eq. m with N = 1 
and initial conditions ip m (t = ^o) = We choose 

{a,p,A} = {10- 8 ,10“ 3 , 6 x 1CT 4 }, so P 2 = 10, and 
Pi = 1.5. In the figure, snapshots of the levels popu¬ 
lation (a,c,e,g) and the electric field (b,d,f,h) are shown 
for multiple r, respectively. The inset also shows the to¬ 
tal energy, '}2 rn 0J rn \il> m \ 2 , which is seen to increase with 
time in a ladder-like manner. The moments of time at 
which m —>■ m +1 transitions occur agree with the theory, 
which predicts r m = mP 2 f5|. The 771 th-level occupation 
numbers, |V’m| 2 , also agree with the theoretical predic¬ 
tions, namely, \4> m \ 2 = TI ^ 1 where Vj^j +1 are 

given by Eq. (USD- For example, at the adopted param¬ 
eters, one has |V> 7 (r = 55 )| 2 = 0.83, which deviates by 
only about 1 % from the value obtained in the simulation. 

A similar argument justifies our neglecting the terms 
of higher order in A in the derivation of Eq. (J5D ■ A di¬ 
rect calculation shows that those terms have the form 
O(c 0 dE 1 A 2 ), so they cause additional “subharmonic” res¬ 
onant transitions m — > m + 2 with the effective driving 
parameter Pi ~ j3A 2 / y/a. The corresponding probabil¬ 
ity P jji — yrn-\- 2 is given by an expression similar to Eq. (USD- 
but Pi should be replaced with Pi 0 - For small enough 
A, one obtains P rn ^ rn+2 ~ (3 2 A 4 /a -C 1. To have both 
fa 1 and P m ^y m+2 <C 1 , one must then require 


FIG. 3: (color online) Langmuir wave evolution in the down- 
chirped AR regime (P 2 = 0.1). The plotted quantities are the 
same as in Fig. [2] but the wave is initialized at larger m, and 
the chirp of the modulation frequency has the opposite sign. 


the condition (which is satisfied in our simulations) 

a 1 ' 2 < A « 0- 1/2 a 1/4 . (11) 

Also note another, kinetic restriction of the LC mech¬ 
anism. It stems from collisionless dissipation, which is 
not contained in our fluid equations. The local Landau 
damping rate for the roth mode in Maxwellian plasma is 
given by y m « k~ 3 exp(-fc“ 2 /2 - 3/2) 0, where 

k m = 77 i/?/3, as before. Therefore, during the transi¬ 
tion between neighboring levels, which occurs on the time 
scale At = (r m+ i —T m )/y/a = /3/a, the energy decreases 
by the factor exp(2r m ), T m = 7 m At. The value of F m 
grows rapidly with m. Hence, just as the wave energy is 
shifted to the mode with m = mkinj defined as that hav¬ 
ing r m ~ 1 , the energy is transferred to electrons almost 
momentarily, heating the tail distribution much like in 
Ref. [25]. At the parameters used in Fig. [l] 77iki n = 9. 
This is larger than the maximum m attained in the sim¬ 
ulation, so neglecting Landau damping is justified. 

AR regime. — In contrast to the LC dynamics, in the 
limit P 2 <C 1, many levels are coupled simultaneously. It 
can be shown then (e^-, using the Wigner phase space 
approach, as in Ref. |5J) that quantum LC continuously 
transforms into the classical AR by decreasing the ef¬ 
fective Planck constant, P 2 . The electric field dynamics 
is then understood as AR acceleration of plasmons that 
satisfy tod = kdV g , where 
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kpt 
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( 12 ) 


Notably, this is the same “group-resonance” condition 
that was recently discussed in Ref. 0 (see Ref. 271 
too). Also notably, the AR acceleration of plasmons 
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that we report here is akin to AR acceleration of res¬ 
onant electrons in phase-mixed nonlinear waves such as 
Bernstein-Greene-Kruskal waves js, ESEi. The AR dy¬ 
namics of Langmuir waves is illustrated in Fig. [2J for 
which we adopted {a,/3,A} = {10 -8 ,10~ 5 ,4.8 x 10 -3 }, 
so P -2 = 0.1. The simultaneous coupling of many lev¬ 
els is clearly seen in the left subplots, while the right 
subplots present the electric field of the autoresonant 
plasmon (solid blue lines) driven by the chirped density 
modulations, n d (red dashed lines). In this case, Landau 
damping affects the total wave energy only by about 1% 
overall, so the kinetic restriction is not essential. 


Interestingly, the dynamics is reversible, at least ap¬ 
proximately. In fact, one can just as well capture a 
wave envelope at large k and then transport it down 
the spectrum, much like a trapped charged particle 
can be decelerated by a resonant field. We demon¬ 
strate the effect in Fig. 0 where the initial condi¬ 
tions are Eg{x) ~ exp(— (x — Xo) 2 /2a 2 )sm(k mr x), 
where m r is the resonant wave number. In this ex¬ 
ample, we apply down-chirped driving phase tpd = 
Um r ,mr+it - at 2 /2 + ip 0 with {/?, x 0 , cr, m r , a, A, <p 0 } = 
{10 -5 ,0,0.05,100,10 -8 ,0.004,7r}, where <p 0 was chosen 
such that the plasmon is initially phase-locked with n d . 

Variational formulation. — Let us now recast our the¬ 
ory in a form that is not restricted to Langmuir waves 
but allows extending the above results to general nondis- 
sipative linear waves. Any such wave, described by some 
real field E(i, x) (electric field being an example), can 
be assigned a Lagrangian bilinear in E, namely, of the 
form |s] 


L = 


E ■ T>(t, x, idt, — *V) • E d 3 x. 


(13) 


The differential operator T> can be considered Hermi- 
tian without loss of generality. Suppose now that T> = 
T> o + T>d, where T> o is some Hermitian operator that 
determines the Lagrangian in stationary homogeneous 
medium, and TQ d is a Hermitian operator that governs the 
mode interaction driven by some weak modulation. Let 
us represent the field as E(i, x) = Re J2 m £m(t) e m (x). 
Here £ m are complex amplitudes and e m are orthonormal 
eigenmodes corresponding to the eigenfrequencies io m of 
T> q. Then, L = m L m + Ld , where the two terms are, 
respectively, due to T>o and T> d - Specifically [3], 

Lm = 2 (VVreVVn l ^m\'4 , m\ : (14) 

where the complex amplitudes are defined such that 
\ifm\ 2 is the action of the mth unperturbed mode; i.e., 
Pm'ifrm = £mi where p m = dZ/. m /dca m , L m = (1/2) J e m • 
T >o ■ e m d 3 x, and w m are found by solving L m (w m ) = 0 


). Also, it is easy to see that 


L d « E 


(15) 


hm,m' — 2 Pm Pm' I ® m * d ' d X. 


Then the equation for ^ m has the form 0 , and, since 
h m ,m' is Hermitian, it manifestly conserves the wave total 
action, For electromagnetic waves in particu¬ 

lar, one has T>d = Xd/^ 71 "); where \ d is the modulation- 
driven perturbation to the medium susceptibility. This 
can be seen, for instance, by comparing Eq. m with 
its geometrical-optics limit [29| • In the case of Langmuir 
waves, x d = J r , and T>q = e/(8ir), where e is the di¬ 
electric permittivity operator; then one recovers Eq. 
Finally, we note that it is straightforward to apply this 
approach to other classical waves [3Cj; hence, our further 
observations of LC and AR apply too. 

Summary. — We report the first theoretical pre¬ 
diction of LC in a classical system. Specifically, we 
show that quasiperiodic chirped modulations of the 
background density can couple discrete eigenmodes of 
bounded plasma to produce controllable shift of the wave 
spectral energy distribution. Apart from academic in¬ 
terest, the new method of continuously controlling the 
wavelength of Langmuir oscillations might find practical 
applications, such as regulating coherent Raman scat¬ 
tering of laser radiation in plasma for generating short 
ultraintense pulses 31]. Our results indicate how similar 
techniques might be practiced with other plasma modes 
too, or in other settings such as waveguides or photonic 
crystals. Our work also bridges a number of effects that 
were previously considered unrelated. In particular, plas¬ 
mon acceleration reported here can be seen as the reso¬ 
nant counterpart of the adiabatic ponderomotive effects 
on waves discussed recently in Ref. [26j. It is also akin 
to the resonant acceleration of charged particles trapped 
by chirped nonlinear plasma waves [16, 28L Finally, our 
results further advance the general idea [1H3] that posing 
classical waves in quantumlike terms can be quite fruitful. 
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